Stability and Ensemble Inequivalence in a Globally Coupled System 
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We consider a system of globally coupled rotors, described by a set of Langevin equations, and 
examine stability of the incoherent phase. The corresponding Fokker-Planck equation, providing a 
unified description of microcanonical and canonical ensembles, bears a few solutions, depending upon 
the ensemble. It is found that the stability of each solution varies differently with the temperature, 
revealing the inequivalence between the two ensembles. This also suggests a physical explanation of 
the quasi-stationarity observed in recent numerical results. 
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A system of oscillators, coupled sinusoidally with each 
other, provides a prototype model for a variety of os- 
cillatory phenomena in nature. With short-range cou- 
pling, i.e., nearest-neighbor interactions, such an oscil- 
lator system describes an array of Josephson junctions 
which has been a subject of extensive studies 0. The op- 
posite limit, corresponding to infinite-range interactions, 
leads to a system of globally coupled rotors, which sim- 
ulates those systems having long-range forces in physics 
and biology. In spite of the mean-field character of such 
a system of globally coupled rotors, the system turned 
out to exhibit rich features in dynamical and statisti- 
cal properties 0j 0, S IE H- Specifically, one can 
find an analytic solution in the canonical ensemble and 
the system has an equilibrium phase transition at a fi- 
nite critical temperature in the ferromagnetic case. On 
the other hand, direct dynamical simulations reveal that 
it takes quite a long time for the system to reach the 
thermodynamic equilibrium [^|^. This anomalous diffu- 
sion phenomenon, dubbed quasi-stationarity, is observed 
only in microcanonical calculation, thus suggesting that 
there may exist inequivalence between microcanonical 
and canonical ensembles. Such a quasi-stationary state 
is believed to survive down well below the equilibrium 
critical temperature, attracting recently much attention 
together with controversy 0, Q . 

This work is to examine the stability of the incoherent 
(paramagnetic) phase in the system of globally coupled 
rotors and to provide an explanation as to the origin 
of the inequivalence between microcanonical and canon- 
ical ensembles. We consider the Fokker-Planck equation 
(FPE), which governs the time evolution of the system 
in both ensembles and thus provides a unified descrip- 
tion, and obtain its stationary solutions, depending on 
the ensemble. The stability of each solution is probed 
and found to behave differently with the temperature. 
In particular, the microcanonical ensemble gives rise to 
additional incoherent solutions, which turn out to be neu- 
trally stable even below the critical temperature. Such 
neutral stability may give rise to quasi-stationarity, thus 
suggesting a resolution of the controversy. 



The dynamics of the coupled rotor system is governed 
by the set of equations of motion for the phase (pi of the 
ith rotor: 



,) = o, 



(1) 



where M is the inertia of each rotor and Jij represents 
the coupling strength between rotors i and j. With the 
introduction of the canonical momentum pi = Mc/ji , the 
above equations are transformed into a set of canonical 
equations with the Hamiltonian 
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on which the microcanonical description is based. 

On the other hand, in the canonical description the 
system is in contact with a heat reservoir of temperature 
T and described, in a most general way, by the set of 
Langevin equations: 



M(t}i + 



i) = Vi 



(3) 



where T is the damping coefficient and the Gaussian 
white noise rjiit) is characterized by the average {r]i(t)) = 
and the correlation {vi{t)Vj{t')) = 2TTS,jS{t-t'). It is 
straightforward from Eq. Q to derive the FPE for the 
probability distribution P{4>i,pi,t): 
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Note that Eq. ^ with P set equal to zero reduces to Eq. 
(0, showing that Eq. Q provides the starting point for 
both descriptions: the microcanonical one (P = 0) and 
the canonical one (P ^ 0). In particular, the stationary 
solution of Eq. Q is given by the canonical distribution 
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P'-°'^{(f>i,p,) cx e"'"/^ with the very Hamiltonian in Eq. 
(0) regardless ofV being zero or not. 

In the case of global coupling = J/N, the set in 
Eq. © decouples into a single-particle equation 



M(t> + T(t> + JA sin((/) ~ 6) ^ r], 



(5) 



where the order parameter Ae' 



N- 



mea- 



sures the emergence of coherence in the system and the 
rotor index i has been suppressed. This in turn leads 
to the FPE for the single-rotor probability distribution 



dP 

dt 
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p + jA sin((/) - 61) 4- rr 



d_ 

dp 



P, 
■(6) 

which reduces, in the absence of damping (F = 0), to the 
FPE for the microcanonical ensemble: 

dP p dP ,^.,,^,dP 

In terms of this probability distribution, the order param- 
eter is simply given by Ae'^ = (e*"^), where (• • • ) denotes 
the average over the distribution P(0,p, i). As pointed 
out for the general case, both Eqs. (jHl and support 
the same stationary solution P^^\(l),p) — Z^^e^^/"^ 
with the single-particle Hamiltonian 



n = 



2M 



JA cos ( 



(8) 



where the overall phase 9 manifesting the U(l) symme- 
try has been absorbed into the definition of </). Note that 
T is given in Eq. ® but remains arbitrary for Eq. Q: 
In the microcanonical ensemble the temperature should 
be defined as a measure of the average kinetic energy ac- 
cording to {p^) = MT. The partition function is deter- 
mined by normalization: Z = {2n)~^ J dp J dcj) e~^/^ = 
V2ttMTIo{x), where /o(a;) is the modified Bessel func- 
tion of the zeroth order with x = JA/T. 

This approach based on the FPE thus provides a uni- 
fied description of microcanonical and canonical ensem- 
bles, and makes clear that both ensembles generate the 
same equilibrium behavior described solely by the 
same stationary distribution p('''((/),p). Namely, in both 
ensembles the equilibrium order parameter is given by 



which, upon expanding e^'^"'"^ in terms of the modified 
Bessel functions, yields 



T 

7" 



h{x) 



(9) 



This determines whether the system exhibits coherence 
(A 7^ 0): The ordered phase emerges when T/ J is smaller 



than the slope of I\{x) / Iq{x) at a; = 0, and the ferromag- 
netic system (J > 0) undergoes a phase transition at the 
critical temperature Tc — J 12. 

We now turn our attention to the stability of the in- 
coherent phase, which depends crucially on the ensemble 
employed, as shown below. For A = 0, the stationary 
solution P^'^'>{cj),p) reduces to the Maxwell distribution, 
describing the incoherent phase. Unlike Eq. ©, how- 
ever, Eq. lO, the FPE for the microcanonical ensem- 
ble, allows (incoherent) solutions of the general form: 
^'"'(0,^) = fo{p), an arbitrary function of p, includ- 
ing the Maxwell distribution. Note that the uniformity 
in (j) guarantees A = since J dcf) e"^ = 0. To probe the 
stability of this stationary solution, we write P{(f),p, t) = 
fo{p) + h{'t>,P.t) and accordingly A{t) = Aq + Ai(t) 
[with Ao = and Ai{t) = J dp J ^ /i((/),p, ^e'^] in Eq. 
(UI), which leads to the stability equation 
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Since fi{(f>,p,t) and Ai{t) are periodic in 
Fourier decompose them into plane waves: 

/i(0,p,t) - E / rft^e'^'=*-"*7fc(P,^) 

7.. 



one can 



Ai(t) 



due 



dpf.i{p,u). (11) 



Note here that the perturbed order parameter is propor- 
tional only to /_i(p, w). Substituting Eq. (|ll|l into Eq. 
(|10ll . one finds the following relations: 



= 



kp 



fk{p,i^) 



M 

dp'f-,{p',uj) 



(4,i-4.-i)f^, (12) 
dp 



satisfied by the Fourier coefficients. Provided lu ^ kp/M, 
it is obvious that all Fourier coefficients vanish except 
those corresponding to k = ±1. When ui = kp/M, we 
have a continuous spectrum which correspond to a time- 
dependent solution of the FPE with A = 0; this will 
be discussed later. For k — ±1, we divide Eq. ((T^ by 
(uj =p p/M) and integrate over p, to obtain 



[1 + X(^)] / rfp/-i(p,c^) = 0. 



where the response function is given by 
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dp 



m 

P + Cj 



(13) 



(14) 



with /o(p) = dfo/dp and u> = Muj. Accordingly, for a 
nontrivial solution, we must have 1 -I- xi'^) — 0- Since 
x(a-') in Eq. (|14|) has a simple pole at p = —Cj on the 
complex p plane on which the integration is performed. 
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one needs to extend the frequency to complex values, 
i.e., (D = (I'r + iCoi- The appropriate analytic continuation 
then leads to the following relation, depending on the 
imaginary part of Cj: 



depending on the sign of uji, leads to 



for oji > 



for tj,- < 0, 



where V stands for the principal part. 

We thus have three different situations to determine 
the frequency: When LOi > 0, the perturbation Ai(t) 
grows indefinitely to make the unperturbed solution 
Ao (= 0) unstable. In this case the conditions for the 
real and the imaginary parts of the frequency to satisfy 
are 
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dp 



{p + OJr)fo{p) 
(p + Ulr)"^ +Ljf 



(15) 



{p + U^r) + UJi 



In the opposite case (coi < 0), the perturbation dies out 
to make the system stable. The condition is then given 
by 
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dp ■ 



{p + UJr)fo{p) 

dp , , - ^2 , -2 + 27rIm/o(-u;) 

(p + Wr ) + ti^j 

' 2^Re/^(-c^) = 0. 



{p + CJr)^ + 

Finally, when = 0, the condition simply reads 

-- 




(16) 



(17) 



and the system is characterized by the neutral stability, 
oscillating with frequency ujr- When this condition is 
met, the system may stay in this state rather a long time 
presumably until higher-order nonlinear terms come into 
play and break neutrality. 

We are now ready to discuss the stability of the in- 
coherent phase associated with various solutions of the 
FPE: the uniform distribution and the rotating distribu- 
tion as well as the Maxwell distribution; the first two exist 
only for the microcanonical ensemble while the third for 
both ensembles. 

Uniform distribution. In this simplest case, mo- 
menta are distributed uniformly in the range [—a, a], 
i.e., /o(p) = l/2a for —a < p < a. Known as a water- 
bag distribution, this has been adopted in most of dy- 
namical calculations performed so far. Substitution of 
/^(p) = (2a)-i[(5(p-Ha) - &{p-a)\ into Eqs. (HHJ) - (113, 



— ^\l -Jr-; ^ f ttI I ujr = Q for a < ai 



= ^V(^) for«>ai, (18) 

which indicates that the frequency of Ai{t) and fi{t) 
has only the imaginary /real part for a smaller/larger 
than ai = yOluT^- Recalling the definition of the 
temperature in the microcanonical description, we have 
T = ju = oP-jiM, from which we find Ti = J/6. At 
low temperatures {T < Ti) the frequency of the perturba- 
tion has the imaginary part, giving rise to the growth (as 
well as decay) of the perturbation. Accordingly the solu- 
tion is unstable. However, at temperatures T > Ti, the 
system oscillates with the (real) frequency, and neither 
dissipation nor growth can be observed. We point out 
here that Ti above which the neutral stability emerges is 
well below the mean field critical temperature Tc = J/2. 
This is precisely what has been observed in various nu- 
merical simulations based on the microcanonical ensem- 
ble: Starting from the incoherent initial configuration 
(Ao = 0) with the uniform distribution, many authors 
reported that the system remains in the incoherent state 
for a long time, before it eventually falls into the coher- 
ent state (Aq ^ 0). Such quasi-stationarity persists down 
well below Tc but the present analysis shows that it 
ceases to exist below Ti. 

Maxwell distribution. We next consider the Maxwell 
distribution, faip) = (27rA/r)-i/2g-pV2A/T^ p^^. > q, 

substituting /q(p) = -{p/MT)fa{p) into Eq. and 
noting f^{—p) = —foip), we find that LUr = is clearly a 
solution of the second one in Eq. 1)15(1 . From the first of 
Eq. p5|l . we obtain 



1 



2T 
T 



y/TTye^ erfc(y). 



(19) 



where erfc(7/) is the complementary error function with 
y = \J Mj^Tuji. Since the right-hand side is monotoni- 
cally increasing from zero to unity as y grows from zero, 
there is one real solution present for T < J/2 and no 
solution for T > J/2. This implies that the system is 
unstable for T < J/2 due to the presence of the solution 
LUi > 0. For < 0, we first note 

from which we see Re/Q((I>r+*Wi) = for J;,. — while 
lnif[^{i^,) = -(27rM3r3)-i/2^^^g<i?/2MT^ rpj^ig suggests 
that (Dr = is again a solution of the second equation in 
Eq. p6|l whereas the first one reduces to 



2T 



1 = V^|y|e^'[2-erfc(|y|: 



(20) 
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Similarly to the Wi > case, the right-hand side of the 
above equation is monotonically increasing with \y\, ren- 
dering a solution uji < only for T > J/2. Finally, when 
uji — 0, we have also uj^ = from /g(— (1)^) = 0; this triv- 
ial solution appears at T = J/2, seen from the first one 
of Eq. ((T7|l . We thus conclude that the Maxwell distri- 
bution, leading to the incoherent phase at T > Tc = J/2, 
becomes unstable below Tc, where coherence develops. 

This conclusion holds also for the canonical ensemble, 
where the Maxwell distribution is the only stationary dis- 
tribution: Equation (jSJ leads to the stability equation in 
the form of Eq. (|10|l with the additional term depending 
on F; this, however, has no effects on the perturbation in 
(j), giving rise to the same instability below Tc. 

Rotating distribution. In addition to the time- 
independent solutions discussed so far, Eq. {Tj) also 
possesses a time-dependent solution of the general form 
P'^°\(l),p,t) = u{(j)-pt/M) for A = 0. In this rotating 
solution the phase grows continuously with a continuous 
spectrum as mentioned previously. Requiring periodicity 
in d), we write 



Fk{p), 



(21) 



where Fk{p) is an arbitrary function except for J±i(p) = 
due to A = 0. The stability analysis for this solution 
is entirely similar to that of stationary solutions and will 
not be repeated here. The response function is given by 



JM [ 



dp 



P 



P 



(22) 



in place of Eq. 1)14(1 . and appropriate changes should 
be made in the subsequent equations through l(T7|l . Al- 
though one can consider various distributions for Fk 0, 
we only quote the results from uniform distributions for 
Fo(j}) and F-2{p), with the distribution width 2a as be- 
fore: 



= ±y(^) uJi^O for a > a2 (23) 

with a2 = J/M. Again, one can relate the average 
kinetic energy with the temperature: T = (jP)/M — 
a'^/3M, from which one has T2 = a^/SM = J/3. Thus 
the rotating solution has neutral stability for T > T2 and 
becomes unstable below T2. Note once again that T2 is 
lower than the equilibrium critical temperature Tc- 

Table I summarizes our results of the stability analy- 
sis for various solutions, depending on the temperature. 
While coherence develops in equilibrium below the tran- 
sition temperature Tc = J/2, the dynamic behavior of 
the incoherent phase depends on the ensemble: In the 
canonical ensemble the Maxwell distribution is the only 





Maxwell 


uniform 


rotating 


J/2 < T 
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n 
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J/3 < T < J/2 
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J/6 < T < J/3 


u 


n 


u 


r < ,7/6 


u 


u 


u 



TABLE I: Stability of the incoherent phase for various dis- 
tributions. Whereas the Maxwell distribution corresponds to 
both ensembles, the other two correspond only to the micro- 
canonical ensemble. The letters s, n, and u stand for stable, 
neutrally stable, and unstable, respectively. 



stationary distribution corresponding to the incoherent 
phase and becomes unstable below Tc. On the other 
hand, the microcanonical ensemble gives rise to addi- 
tional distributions such as the uniform and the rotat- 
ing one, which remain neutrally stable even below Tc- 
Putting these together, we are led to the plausible conclu- 
sion that the quasi-stationarity found in numerical sim- 
ulations below Tc has its origin in the neutrally stable 
solutions in the microcanonical ensemble and will not 
appear in the canonical ensemble. 

Similar calculations for the antiferromagnetic system 
(with J < 0), where coherence does not emerge [see Eq. 
©], lead to that all three types of solutions are always 
neutrally stable 9]. This may give a clue to nature of 
the non-stationary behavior observed recently |lCll| . 
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